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Abstract. We investigate the plausibility of some models emerging from an 
algorithm devised to generate a one-parameter family of interior solutions for the 
Einstein equations. It is explored how their physical variables change as the family- 
parameter varies. The models studied correspond to anisotropic spherical matter 
configurations having a non local equation of state. This particular type of equation 
of state with no causality problems provides, at a given point, the radial pressure not 
only as a function of the density but as a functional of the enclosed matter distribution. 
We have found that there are several model-independent tendencies as the parameter 
increases: the equation of state tends to be stiffer and the total mass becomes half 
of its external radius. Profiting from the concept of cracking of materials in General 
Relativity, we obtain that those models become more stable as the family parameter 
increases. 
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1. Introduction 

Compact objects are one of the most fascinating entities known in our Universe. These 
very particular astrophysical bodies (white dwarfs, neutron stars, quark stars, hyperon 
stars, hybrid stars or/and magnetars) seems to be rehcs of luminous stars which, in turn, 
are considered to be born from extreme events such as in supernova explosions. These 
objects are thought to have cores with densities higher than nuclear (p ~ (10 — 20) po 
with Po ~ 2 X 10^^ grcm~^) and enormous gravitational binding energy ~ 5 x 10^^ 

ergs ~ 0.2Mc^). 

The true equation of state that describes the properties of matter at such a 
high densities is essentially unknown [l]-[2]. Currently, what is considered to be best 
understood in this active field, comes from the experimental insight and extrapolations 
from the ultra high energy accelerators and cosmic physics (see [3] and references 
therein). Having this uncertainty in mind, it seems reasonable to explore what is allowed 
by the laws of physics within the framework of the theory of General Relativity. 

We shall consider a general relativistic matter spherical matter configuration, 
having a particularly Nonlocal Equation of State [NLES from now on): 

P(r) = p(r) - ^7 r f^pif) df ^ pir) - 3P(r) + r [p(r) - P(r)]' = 0,(1) 
r-^ Jo 

which is the static limit of a more general relation between two components of the energy 
momentum tensor: 

-T} = T°-- /'T°df ^ [t° + 3Tj + r fT° + Ti' 

r Jo L ^ ^ 



0,(2) 



where u is the proper velocity in the radial direction and the prime denote derivative 
with respect to the radial coordinate r (see [4, 5] for details). 

Equation (2) was obtained when constant compactness (i.e. constant gravitational 
potentials at the surface of the matter configuration) was required for imperfect 
anisotropic fluid (unequal stresses, i.e. P P± and heat flux) matter configurations. It 
was also found that these type of anisotropic configurations admit a Conformal Killing 
Vector Field (see reference [4] for details). 

It is clear that in equation (1) a collective behavior on the physical variables p(r) 
and P{r) is present. The radial pressure P(r) is not only a function of the energy density, 
p(r), at that point but also its functional throughout the rest of the configuration. The 
inner distribution of matters to a given point, and not only the density at that point, 
contributes to the value of the radial pressure. It is worth mentioning that there is no 
possible causality drawback with equation (1) because it is the static limit of Equation 
(2). When the evolution of a density perturbation is examined, the full dynamical 
Equation (2) has to be considered. 

This type of equation of state, originally proposed by D. G. Ravenhall and C. J. 
Pethick in 1994 [6], has proven to be very fruitful describing a variety of relativistic 
astrophysical scenarios [7, 8, 9, 10]. Particularly, in [7] the authors, considering 
the framework of the non commutative geometry, describe a mini black hole having 
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a Schwarzschild geometry outside the event horizon, but whose standard central 
singularity is replaced by a self-gravitating droplet. The energy- momentum tensor of 
this droplet is an anisotropic fluid obeying a NLES, which allows to have a droplet with 
a finite radius and positive definite pressure distribution at the interior. More recently, 
the concept of quasi-local variables was introduced which leads to a simplified criterion 
of stability of matter configurations based on the profiles mass M vs radius R [9]. 

In General Relativity, static spherically symmetric perfect fluid distributions are 
described by three independent Einstein equations for four variables (two metric 
functions, the energy density and the isotropic pressure). As in the case of classical 
hydrodynamics, in order to integrate the system, additional information has to be 
provided in the form of an equation of state or an heuristic assumption involving metric 
and/or physical variables. This situation suggests the possibility to obtain any possible 
solution by giving a single generating function, inspired in some physical intuitions. A 
method to procure solutions in this way has been recently presented by K. Lake [11] 
for local spherical isotropic fluid (equal stressed along radial and tangential directions). 
This algorithm was extended by Herrera and collaborators [12] in order to consider 
anisotropic spherical distributions of matter. 

In the literature there are several proposed strategies (with a variety of heuristics) 
to produce families of solutions having a parameter that generates new family-members 
(see, for example [13]). It is interesting to inquire about how the physical variables 
change as the parameter is altered and, more specifically, how the variation of the 
family-parameter affects the stability of the corresponding matter configuration. If the 
algorithm to generate new solutions leads to plausible models of compact objects and, 
as the family-parameter varies the family members become better models, thus the 
proposed method will have more physical meaning. In this work we shall investigate 
families of NLES anisotropic matter configuration generated with the Lake-Herrera 
method by exploring the plausibility of their physical variables and by studying their 
potential stability when the parameter is changed. 

The structure of the present paper is the following. The next section sketches 
the method to generate families of exact solutions having NLES. In Section 3 we shall 
workout three examples of families of matter configuration having NLES. We close the 
work with a brief discussion of the results in the Section 4. 

2. The Method and the NLES 

In order to explore the feasibility of NLES in the context of the Lake-Herrera method, 
we shall use an energy-momentum tensor represented by 



to describe static, spherically symmetric, anisotropic (unequal stresses, i.e. P ^ P±) 
bounded configurations in General Relativity. 



Following [14], we can build a set of four independent scalars as: G^M°n^, 
Gl = G^u'^up, G2 = G^n'^np and G = G^. Where G^ is the Einstein tensor, with 




(3) 
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a congruence of unit timelike vectors = (m^,m^,0, 0) and its unit normal field 
satisfying = 0, and inf^np = 1. 

Thus, with the above structure of the energy-momentum tensor, we have the 
following equations relating those scalars and the physical variables of the configuration: 
G = 8tt{p-P- 2Px) , Gl = Svrp , G2 = 8ttP , -G + Gl - G2 = IGirP^ . 

For the perfect fluid case, (P = Pi), we can write an equation involving only the 
above mentioned scalars: 

- G + Gl - 3G2 = . (4) 

Considering an static spherically symmetric space time, described by the line 
element 

ds^ = e'"^"-^ dt' - e2"(^) dr' - (dO^ + sen^^ d<p') . (5) 

If e"^'^*^''^ = 1 — 2m{r)/r then, equation (4) can be written as 

(r — 2m) u" — [u'Y — ru [r — 3m + rm] — rm! + 3m = . (6) 

It is clear that for a given i^(r), the equation can be integrated to obtain m(r), making 
it possible to build an algorithm to generate families of solutions of Einstein equations 
in the case of perfect fluid [11]. The function v{r) is not completely arbitrary because it 
must fulfill a series of requirements that prevent singularities on the invariants obtained 
from the Riemann tensor. 

For a perfect, static and spherically symmetric fluid, the fact that the central density 
pc and the central pressure Pc are finite ensures that all the invariants of the Riemann 
tensor, are regular at the center of symmetry. 

We wish to consider the feasibility of using a NLES to obtain families of exact 
solutions of the Einstein equations through the above mentioned algorithm and to 
explore the behavior of the physical variables as we vary the parameter that generates 
the members of the family. 

The metric corresponding to a static spherical matter configuration having a NLES 
can be written as 

ds^ = e^^^'^e^^ df - e^^^'^ dr^ - [dO^ + sen^^ dcj>^) , (7) 

where A(r) = z/(r) — k and k is a constant [4, 5]. This metric recalls the so called 
isothermal coordinates system [15] which, in turn, is a particular case of the more 
general "warped space-time" (we refer the reader to [16] and references cited therein for 
a general discussion). 

With the metric (7) and the energy momentum tensor (3), the Einstein equations 
can be written as 

1 - e"2^ (1 - 2rA') 
1 - 

A" 
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Now, combining (8) and (9) it results that 

m(r) = 2Tir^{p-P) . (11) 

From the matching conditions, at the boundary of the distribution r = R, we get: 

m{R) 



K 



—2 Ar = ln(l — 2/i), with /i 



R 



M 

'r 



(12) 



Clearly, from equation (1) evaluated at the center of the matter distribution, it follows 
that pc = 3Pc- Additionally, in order to have P{R) = we must satisfy: 



2i?A', 



R 



where \r = X{R) and 



^2Xn ^ 

__ dA(rO 
dr 



(13) 



r=R 



2.1. A family of matter configurations with NLES 

Restating a similar algorithm as the one presented in references [11, 12], we can generate 
families of exact solutions, for the case of NLES, by providing a metric function z/(r). 
Our version for the NLES Lake-Herrera method involves the following steps: 

(i) Provide a metric function i/(r) of the form 

i/(r) = A^$(r) , A^>1, (14) 

where is an integer parameter that defines the family and $(r) is a function that 
must be monotonically increasing with a regular minimum r = and have to fulfill 
with the matching conditions at the boundary of the distribution [11]. 

(ii) From the condition (13) it follows that 



,2 k 



;i - 2^iy 



2 RN<^'ji ' 



(15) 



with = $(i?) and 



d$(r) 
dr 



r=R 



(iii) From the Einstein equations (8)-(10) we obtain the physical variables: 



P 



P 



2rN^' 



2 RN^', 



P^ = — 



-1 + 



R 

2riV$' 



87r 1 + 2 RN(^' 



1 + 2 RN^'j^ 

„2JV(*fl-<I.) 



^2Af($fl-$) 



(16) 
(17) 



(iv) The mass function, equation (11), is: 



32Ar{$H-(l>) 



1 



2 RN^'r 



R 



(19) 



(v) We can also calculate algebraic expressions for the velocities of sound, both radial 
an tangential as 



dP 
dp 



^^2 _ 2r$' - 


1 ' 








vl/r2 + i 




-2R 



(20) 



6 



dp 2 ^r2 + i, 



N 



TV 



where ^ = $" - 2iV$'2. 



2.2. How this family of NLES behaves for Large N ? 

The above method allows us to explore the behavior of the physical variables for different 
values of A^, in particular we can see that for extreme situations, for large values of A^, 
it is clear that: 

Ar»i ^p^Z|J_e2^(*«-*) and P ^ ^e^^^^^"*) , (22) 

R^R 

thus, the material become very stiff, i.e. when ^ 1 the equation of state tends to be 
p = P, since we know that $r — $ > 0. 

Considering the anisotropy of the material 

A^ > 1 — ^ ^ . 23 

p 2r$' ^ ' 

As it emerges from the hydrostatic equilibrium equation for general relativistic 
anisotropic fluids, the gradient of pressures is steeper for positive values of the anisotropy 
(Pj_ — P). Thus, if ($"(r) — 2r$'(r)) > the pressure towards the center will increase 
more rapidly. It is clear that the effect of the change of the anisotropy as A^ increases 
depends on the sign of the equation (<I>"(r) — 2r$'(r)) which could be translated into 
restrictions on the seed function $(r). This will be clear with one of the examples 
worked out in Section 3. 

Now, by evaluating the equation (19) at the boundary r = R, we obtain for the 
total mass M: 

A^ > 1 ^ M ~ ^ , (24) 

regardless of the seed function $(r). 

Finally, the sound velocities when A^ ^ 1 tend to behave like 

2 ($"-2A^<l>'2)r2-2r$' ^ $'"_2A^<I>"$' 

($"_2A^$'2)r2 ^ ($//_2A^$'2)^2 • ^ > 

Which are model (seed function) depend. Next section will be devote to present some 
examples of the above algorithm. 

3. The NLES families by some examples 



In this section we shall workout three examples of families of matter configuration that, 
having NLES, emerge from the method above mentioned. 



7 



3.1. Tolman IV families 



The Tolman IV isotropic static solution was originally presented by R.C. Tolman in 
1939 [17] and it is also found as a particular case of a more general family of solution 
in [18] and [19]. It is, in some aspects, similar to the equation of state for a Fermi gas 
in cases of intermediate central densities. 
Let us consider the a seed function as: 



- In a + — 



(26) 



where a and b are constants. 

From (8) - (10) we get the corresponing physical variables: 



P 



{ab + r^) 



{ab + r 



1 



{ab + r^) 



2^{N+l) 



2k 



6^e2«(2A^ + l)r2 



87rr^ 
Nb^ {ab - r2) 



{ab + r 



2n(7V+1) 



87r (a6 + r2) (^+2) • 
The mass function can also be written as 



m{r) 



(27) 
(28) 
(29) 

(30) 



{ab + r^) 

Now, the two constants are determined from regularity condition at the origin r = 
for the density and pressure, p and P, i.e. 



,2k \ N 



and the constant b is obtained by solving (15), i.e. 



N 2k 



ab 



2\{N+l) 



(31) 
(32) 



i?^(2A^ + 1)] = (ab + R^ 

Another way to see the regularity at the center of these models is to calculate 
central density which, in this case can be written as 



Pc 



2(JV+1) 

Sirb e N ' 



(33) 



substituting the corresponding values for a and b in the expression for m{R) = M 
Equation (30), we obtain 

/ 37V \^ 

1 - 



2 



(34) 



\3N + 8ttR^pJ 

Below, in Table 3.1, there are several cases for different values of A^, where the 
constants a and b are obtained from Equations (31) and (32). For > 4 the calculations 
will depend on the ability to solve analytically equation (32) but, in any case, it is always 
possible to find its numerical solution. 
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It is important to note that that = 1 corresponds to a solution Tolman IV-hke, 
which is anisotropic and has a NLES. The resulting equations of state for values of 
> 1 have completely different properties, and can appreciated from Figure 1. 
Table 3.1 shows some numerical values for compact objects with radius i? = 10 km. 



Model 


/i 


M(M0) 


a 


b/M 


PR X 10'' (t) 


Pc xl0i5(t) 


N = 1 


0.25 


1.69 


0.25 


256 


5.36 


1.61 


N = 2 


0.29 


2.00 


0.41 


122 


6.32 


1.81 


N = 3 


0.31 


2.12 


0.52 


97 


6.72 


1.87 


N = A 


0.32 


2.19 


0.59 


91 


6.94 


1.91 


Ta 


ole 3.1 



The equation of state P = P{p) for some values of are displayed in Figure la, 
with the original Tolman IV like equation of state for the case A^ = 1. We can appreciate 
a satisfactory behavior of the three profiles for the physical variables in Figure lb and 
also form the anisotropy (Pi — P)/p for different values of the parameter A^. 

3.2. Families generated from conformally flat solutions 

Now we shall present two families of models generated from anisotropic conformally flat 
seed solutions. The two conformally fiat seed solutions considered in this section were 
presented by B.W. Stewart [20]. In order to build the family of nonlocal solution, we 
shall proceed in two ways. First (Case 1), we will select as a seed function, one inspired 
by a static conformally fiat solution and we turn it to be nonlocal by the use of the 
method. Secondly (Case 2), we will choose a static solution which is both, conformally 
fiat and nonlocal. In both cases only for A^ = 1 the configurations are conformally flat 
and, as it will become evident, when A^ increases the behavior of the physical variables 
for these to models are qualitatively different. 
Now, considering [20] we have 

$'(r) = - (r^ - 2rm{r))"'^^ , (35) 

ClS Sb condition to be satisfy by the metric function 

eM0=e2*(^ + 5^2^-2*) 2^ (36) 

and the quantities a and b are constants to be determined by boundary conditions. 

Following [20], we start with some known static metric function, m(r), and by 
integrating equation (35) the other metric function is obtained. This is the most critical 
part of this approach and from the four examples shown in [20], two of them have to do 
with solutions by series for $(r). 

f In gr.cm^'^ 



3.2.1. Case 1: Consider Example 1 of [20] and take mass function shown there 



m[r] 



r 
2 



1 - 



A = In 



sin(i^r) ) ' 

where i^' is a constant. From (35) it is obtained 

e*(^) = ^i^rcot {^)+C, 

and if e2*(o) = 1 then C = 0. 

For simphcity, here we will make the following change of variables 



(37) 



(3^ 



/(r) = cot ( ^ 



/' = -f (1 + /^ 



(39) 



We want then to study the feasibility of the method taking the following family of 
metrics 



In 



Klifr 



with n = a + 



4& 



and a = e 



k/N 



(40) 



The field equations arising from the nonlocal condition are 



P 



N [Klip (2/ - Kr{l + p)) + 8r6 (1 + p)] 



and 



KUp 
Stt '' 



P = p~ 

and the mass function 
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m{r) = - 

From the boundary condition P{R) = is possible to obtain the constant b: 







27V" 


1 - A^e^^ 


1 






KUrf 





NKR (1 + F^) - (2iV + 1 - e"'') F 
NKR (1 + F2) + (2iV + 1 - e-«) F 



where F = cot f-^)- From m{R) = M it is obtained 



V 2 y 
ARb = KF 



2 - RRFe""'^^ 



(41) 



(42) 



(43) 



From equations (42) and (43) there results the following transcendental equation for the 
factor KR 



2KRFe^ = ± 



\ 



(4 [2N + 1- e-'^] +N)+N 



1 . 



(44) 



N[l + F2] 

For a fixed value of KR and different values of A^, it is possible to obtain a numerical 
solution of this equation for the amount /i. The constant b is then obtained for (42). 
Following, the central density can be calculated for all 

N 



Pc 



levr 



Ube-^/'' - K' 



(45) 
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By combining this equation with (41) evaluated at the surface is: 

/ 6NKF \ 

1 - 



2 



R[NK^ (3F2 



167rpc 



(46) 



The following table shows some numerical values for a compact object of radius 
= 10 km and KR = 0.1 



Model 


fi 


M(M0) 


a 




PR X lOi^(t) 


Pc xl0i5(t) 


N = 1 


0.30 


2.00 


0.41 


0.02 


6.33 


1.81 


N = 2 


0.32 


2.19 


0.59 


0.02 


6.94 


1.91 


N = 3 


0.33 


2.27 


0.69 


0.02 


7.17 


1.95 


N = A 


0.34 


2.30 


0.75 


0.01 


7.29 


1.97 



Table 3.2 

Figure 2 shows a set of graphs for models here referred to as the Case 1 of 
conformally flat nonlocal configuration. The equation of state P = P{p) for some values 
of N can be appreciated in Figure 2a, as well as the equation of state P = p when N 
= 1 which is the original member of the family. We can also appreciate (Pi — P)/p for 
each N in Figure 2b. 



3.2.2. Case 2: A completely different treatment of the above is to start from 

n 2 



.2A 



2^2 



1 + cot^ {Kr + A) 



Kr 



sin {Kr + A) 



(47) 



where K and A are integration constants. 

The constant A should vanish because it is necessary that the metrical elements 
are regular functions at the origin, i.e. 

A = ^ lim e^^ = 1 . (48) 

It is clear that, by imposing another condition on the metric elements, it is possible to 
completely determine the system of equations. 
Now reapplying the algorithm we get 



$(r) 



In 



sm.{Kr) 

Again, in order to e^'^*-'^-' = 1, we must have A(r) = — Nk . 
From the field equations we get 

NKr sin (2 J-Cr)" 



(49) 



P 



1 

N 



2N-1- 
2N+1 



sm^ (Kr) 
NKr sm{2Kr) 



sin^(-ft'r) 



- 1 



sin^(Kr) 

sin(i^r) 
Kr 




2N 



2N 



2N 



(50) 
(51) 
(52) 



f In gr.cm 
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and the corresponding mass function 



m{r) 



' sm{Kr) 
Kr 



2N 



At the surface where m{R) = M, we should satisfy the following identity 



'sm{KR)\ 



2N 



KR 



sm{KR) = KRe^ 



and the other condition at the surface, P{R) = 0, is 



2N 



NKR sm{2KR) 
sin^ (KR) 

But by using (54) and (55) it is obtained 



smjKR) 
KR 



2N 



K^R^ 



2iV + 1 - e"^ 
2N 



(53) 



(54) 



(55) 



(56) 



This expression must be substituted into the transcendental equation (54) in order to 
obtain (numerically) a value of the mass-radius ratio fi for different values of N. 
Then, the central density turns out to be 

NK^ , , 

Combining this last equation with (53), evaluated at the surface, it results in 



2 



N 



I-kR^P^ 



sm 



N 



R 



cos 



N 



R 



2N 



{51 



The following table shows some numerical values for a compact object of radius 
R=10 km. 



Model 




Me 


KM 


PR X 10^^(1) 


Pc X 1015(1) 


N = 1 


0.396 


2.683 


0.79 


8.49 


2.14 


N = 2 


0.374 


2.535 


0.52 


8.03 


2.07 


N = 3 


0.368 


2.494 


0.42 


7.90 


2.05 


N = A 


0.366 


2.475 


0.36 


7.84 


2.04 



Table 3.3 

It is clear that for values of 7^ 1 Figure 3 shows a set of graphs for models here 
referred to as Case 2. The equation of state P = P{p) for some values of are sketched 
in Figure 3a, as well as the equation of state P = p. We can appreciate the anisotropy 
(Pj_ — P)/p for different A^ in Figure 3b. 



f In gr.cm 



12 



3. 3. The stability of the models 

Fluctuations in density and anisotropy induce total radial forces which, depending on 
their spatial distribution, may lead to the cracking, i.e. radial force directed inwards, or, 
overturning, directed outwards, of the source [21]. We can evaluate potentially unstable 
regions within anisotropic models based on the difference of the propagation of sound 
within the matter configuration. Those regions where the radial sound speed is less 
than the tangential sound speed: v1 < f could be potentially unstable, but if > v'^j_ 
everywhere within a matter distribution, no cracking will occur and the configuration 
will be potentially stable [22]. 

By means of the equations (20) and (21) is possible to determine the sign of 
Sv"^ = — v"^ and to explore how it changes as varies. Figures Ic, 2c, and 3c show 
different profiles of 6Vg for Tolman IV like. Case 1 and Case 2 models, respectively. 

Notice that the original seed functions {N = 1) for models Tolman IV like and Case 
1 are potentially unstable, this means that there are regions where cracking could occur, 
however, as increases the possibility of occurrence of cracking vanishes. Unlike this 
situation the Case 2 is potentially stable for all A^. 

4. Conclusions 

In this work we have found three new families of solutions describing nonlocal anisotropic 
compact objects by applying the Lake-Herrera [11, 12] algorithm. As it was illustrated 
with the physical quantities shown in the above tables 3.1, 3.2 and 3.3, 

0.25 < fi< 0.39 I.6OM0 < M < 2.68M0 1.61 < Pc x 10^^ < 2.04 

the models worked out describe plausible compacts relativistic objects. 
With the particular ansatz (14) considered, we have found that: 

• if A^ increases the material become stiffer and the equation of state approaches the 
limit p = P. This is clear from (22) and also it is apparent from the Figures la, 2a 
and 3a. 

• the effect of the change in the anisotropy is model dependent. It depends on the 
sign of the factor ($"(r) — 2r$'(r)) which could be translated into restrictions on the 
seed function $(r). This model dependent behavior could be appreciated comparing 
figures 2b and 3b. In the first N increases the anisotropy (Pi — P)/p, for 
a particular point within the matter configuration, increases but, for the Case 2, as 
A^ increases the anisotropy decreases. 

• concerning the total mass, M = m{R), it is clear that for A^ ^ 1 we got M ~ ^ 
regardless of form of the seed function $(r). This is consistent with the graphs 
sketched in figures Id, 2d and 3d for the three equations of state considered. It is 
found that total masses approaches a limit lower than 3.4 solar masses, i.e. R/2 V 
N. 
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• if the initial model of the family presents potential instability, it could become 
potentially stable as N increases. This is clear from figures Ic and 2c where the 
relation Svg = — for different values of are displayed. As it can be 
appreciated form Figure Ic 6v^ < for > 2 and from Figure 2c Sv^ < for 
N >1. 

For the considered ansatz (14), the Lake-Herrera method seams to be a very 
reasonable algorithm to generate models for compact object with a NLES. 

It is worth mentioning that the set of nonlocal solutions obtained with this 
algorithm describes physically reasonable compact objects. There is no causality 
drawback for any of the NLES considered, and in fact there is no causality problems for 
any NLES. The Equation (1) is the static limit of Equation (2) and when the evolution 
of a density perturbation is examined, the full dynamical Equation (2) has to be taken 
into account. 
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Figure 1. Tolman IV like 




Figure 2. Case 1 



